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Afcttrac* 

[This  report  analyses  the  performance  of  the  nonlinear  altitude  tracker  developed 
for  the  Active  Beacon  Collision  Avoidance  System  (BCAS)  by  Lincoln  Laboratory. 
[The  tracker  is  intended  for  use  in  the  collision  avoidance  logic  of  the  BCAS 
system.  The  nonlinear  Cracker  performance  was  characterized  through  comparisons 
with  the  previous  altitude  cracker.  The  original  tracker,  used  by  the  collision 
avoidance  logic,  was  a  simplistic  Alpha-Beta  (m-^-0)  tracker. 

The  nonlinear  tracker  performance  evaluation  was  conducted  in  three  phases:  (1)  The 
stand-alone  error  characteristics  of  the  tracker  were  obtained.  Simulated  mode  C 
report  sequences  were  provided  directly  to  the  tracker.  (2)  The  nonlinear  tracker 
was  integrated  directly  into  the  collision  avoidance  logic.  With  the  use  of  the 
Fast-Time  Encounter  Generator  (FTEG),  a  comparative  study  of  performance  with  the 
nonlinear  tracker  versus  the  ^  -  C  tracker  was  made.  (3)  Selected  live  flight 
test  encounters  were  used  to  analyze  the  relative  performance  of  the  a  -  0  tracker 

versus  the  nonlinear  tracker.  •.  n 

fhr  .  -  .*> 

The  stand-alone  analysis  revealed  that  the  nonlinear  tracker  consistently  had 
smaller  maximal  errors  ’in  vertical  race  estimation  and  a  smaller  transient  rate 
response  delay  than  did  the  a  -  0  tracker.  Both  the  live  flight  test  encounter 
simulations  and  the  FTEG  scenario  simulation  indicated  that  nonlinear  tracking 
often  caused  an  increase  in  separation  for  encounters  with  vertically  accelerating 
threats  and  reduced  occurrences  of  inccrrect  command  sense  choice. 


1?.  Kay  W«i|, 


Nonlinear  Altitude  Tracker 
Alpha-Beta  Tracker 

Fast-Time  Encounter  Generator  (FTEG) 
Beacon  Collision  Avoidance  System  (BCAS) 


IS.  OlttflkvNan  Imurnil 

Document  is  available  to  the  U.S.  public 
through  the  National  Technical  Information 
Service,  Springfield,  Virginia  22161 


Farm  DOT  F  1700.7  <a-7j> 


Security  Claitif.  (•(  *Ki»  papa) 


•FRY  T  1 1I4I-? 


Rapraductian  af  camplatad  pofa  sutKaritaJ 


22.  Prica 


ACKNOWLEDGEMENTS 


The  author  wishes  to  express  appreciation  to  Richard  Cleary  of  the  Federal  Aviation 
Administration  Technical  Center.  Mr.  Cleary  was  instrumental  in  identifying  flight 
test  Active  Beacon  Collision  Avoidance  System  tracking  anomalies  discussed  in  this 
report.  The  author  also  wishes  to  express  appreciation  to  John  Andrews  of  Lincoln 
Laboratory,  the  deve'.oper  of  the  nonlinear  tracker  analysed  in  this  report. 
Mr.  Andrews  showed  extreme  patience  in  listening  to  and  answering  questions  about 
the  tracker. 


Aooesalon  For 

FT IS  OKA* I 

DTIC  TAB 

X 

□ 

Unannounced 

Just If teat Ion _ 

□ 

L  _ 

_ 1 

I  Distribution/ 

I  Avallallllty  Codes 

Avail  and/or 

Dlat 

Jli 

Special 

DTIC 

SELECTED 
DEC  28  1981  ,  ■ 

D 


TABLE  OF  CONTENTS 


Pag* 

INTRODUCTION  1 

Purpose  1 

Scope  1 

Background  1 

DISCUSSION  2 

Previous  Vertical  Tracker  2 

Initial  Flight  Test  Results  3 

New  Nonlinear  Tracking  Approach  6 

RESULTS  AND  ANALYSIS  8 

Stand-Alone  Steady  State  Performance  8 

Acceleration  Performance  10 

Deceleration  to  Level  Flight  Conditions  19 

Deceleration  to  Other  than  Level  Flight  Conditions  22 

Pitchover  Maneuvers  22 

Stairstep  Maneuver  27 

Fast-Time  Encounter  Generator  29 

Review  of  Live  Flight  Test  31 

CONCLUSIONS  36 

RECOMMENDATIONS  38 

REFERENCES  39 

APPENDIX 


v 


LIST  OP  ILLUSTRATIONS 


Figure  Pag* 

!  Rat*  ard  Position  35-Second  Projection  Error  Cycles  4 

(I  -  0.10) 

2  Level  Plight  Mod*  C  Report  Excursion  Patterns  6 

3  Mode  C  Excursion-Induced  Tracked  Rat*  Errors  7 

4  Comparison  of  Limit  Cycle  Oscillation  for  500  Pt/Min  Climb  9 

5  Maximal  Error  Comparison  for  Pinal  Steady  Star*  Climb  10 

6  Comparison  of  Acceleration  Performance  —  Acceleration  12 

Race  a  0.25  g,  Final  Rat*  •  1,000  Ft/Min 

7  Comparison  of  Acceleration  Performance  —  Acceleration  13 

Rat*  a  0.25  g,  Final  Rate  a  2,000  Ft/Min 

8  Comparison  of  Acceleration  Performance  —  Acceleration  14 

Rate  a  0.25  g,  Final  Rate  a  3 >000  Ft/Kin 

9  Slow  Acceleration  and  Nonconataat  Bin  Occupancy  Performance  13 

Comparison 

10  High  Acceleration  and  Nonconstant  Bin  Occupancy  Performance  17 

Comparison 

11  Deceleration  to  Level  Flight  Condition  Performance  g  a  -0.1  20 

12  Deceleration  to  Level  Flight  Condition  Performance  g  a  -0.5  21 

13  Comparieon  During  Deceleration  to  Fixed  Rate  Other  Than  24 

Level  Flight 

14  Pitchover  Maneuver  Perfonunce  Comparison  25 

15  Typical  Stairstep  Maneuver  27 

14  Stairstep  Maneuver  Performance  Comparison  28 

17  Unequipped  Intruder  Encounter  Sequence  Which  Results  in  30 

Incorrect  Sense  Choice 

18  Geoauttries  Used  in  Comparison  of  Detection  Performance  31 

19  Vertical  Profile  of  DC-9  Encounter  34 


LIST  OF  ILLUSTRATIONS  (Continued) 


Figure  Pege 

20  Live  Flight  Test  Encountei  —  Alpha-Beta  Tracking  Results  35 

21  Live  Flight  Teat  Encounter  —  Nonlinear  Tracking  Reaulta  35 

22  Tracker  Reaponaea  for  DC-9  Vertical  Profile  36 

LIST  OF  TABLES 

Table  Page 

1  Quantisation-Induced  Cyclic  Peak  Error  in  the  Rate  Eatiiaate  5 

and  the  Impact  on  the  Projected  Poaition  g  ■  0.10 

2  Trackera  Response  Delay  as  a  Function  of  I  7 

3  Comparison  of  Time-Dependent  Average  Error  Magnitudes  During  IS 

Accelerations  —  Ej  in  Ft/Sec 

A  Comparison  of  Time-Dependent  Average  Error  Magnitudes  During  23 

Decelerations  to  Level  Flight 

5  Comparison  of  Tracker  Response  Delays  Following  -0.5  g  26 

Pitchover  Maneuver  to  -1,000  Ft/Min 

6  Reduction  in  Duration  of  Incorrect  Sense  Choice  Periods  in  32 

Seconds  Achieved  with  Nonlinear  Tracking 
(-0.25  g  Deceleration) 

/  Reduction  in  Duration  of  Incorrect  Sense  Choice  Periods  in  33 

Seconds  Achieved  with  Nonlinear  Tracking 
(-0.5  g  Deceleration) 


INTRODUCTION 


PURPOSE. 

The  purpose  of  this  report  ia  to 
describe  the  analysis  of  the  nonlinear 
tracker  algorithaa  developed  by  Lincoln 
Laboratory  (reference  1).  The  nonlinear 
tracker  was  developed  to  replace  the 
current  alpha-beta  (  a  -  I  )  tracker  in 
the  Active  Beacon  Collision  Avoidance 
Syatea  (BCAS)  intruder  and  own  aircraft 
tracking  logic  (reference  2).  Tne 
report  documents  the  performance  of  the 
nonlinear  tracker  in  teras  of  vertical 
rate  error  characteristic  and  vertical 
separation  at  the  closeat  point  of 
approach  (CPA)  following  BCAS  coeaunds. 

SCOPE. 

This  report  is  intended  to  characterise 
the  iapact  of  the  mode  C  altitude  report 
quantitation  on  vertical  rate  tracker 
estimates.  This  report  presents  results 
of  last-tiae  simulations  of  nonlinear 
tracker  perforaance. 

An  analysis  was  stade  of  the  stand-alone 
vertical  rate  and  position  error  charac¬ 
teristics  and  coapariaon  with  the 
perforaance  of  the  •-  g  vertical 
rracker  using  specially  developed 
conflict  scenarios  and  encounters  based 
on  live  BCAS  flight  test  data. 

BACKGROUND. 

The  BCAS  logic  requires  the  accurate 
sequential  estimation  of  own  aircraft's 
and  the  intruder's  vertical  position  ard 
vertical  race.  This  information  is 
generated  every  logic  cycle  (nominally 
every  second)  in  the  own  airersfe  and 
intruder  aircraft  tracking  modules.  In 
Che  preaence  of  Biasing  target  reports, 
the  tracking  modules  can  coast  the 
vertical  poaition  and  vertical  race 
based  on  previous  track  history.  On 
every  logic  cycle  after  the  logic  has 
updated  own  and  intruder  aircraft 
tracks,  the  vertical  poaition  and  race 


estimates  are  uaed  to  perfora  up  to 
three  vital  collision  resolution 
functions. 

The  first  function  is  to  determine  which 
intruders  are  threats  and  iJ  a  reso¬ 
lution  advisory  is  required.  Previous 
testing  of  the  BCAS  logic  showed  that 
the  errors  associated  with  the  a  -  I 
tracked  vertical  poaition  and  vertical 
rate  only  reused  a  1-  to  5-second 
variation  in  the  tiae  at  which  the  BCAS 
coanand  first  appeared  (i  ference  3). 

The  second  use  is  to  determine  the  sense 
of  escape  aaneuver.  Active  BCAS  can 
only  provide  escape  resolution  in  the 
vertical  doaain.  The  sense  of  aaneuver 
ia  selectad  by  modeling  an  own  aircraft 
response  to,  first,  a  cliab  coaaand 
and,  then,  a  descent  coaaand.  The  sense 
of  xar.euver  which  would  provide  the 
greater  separation  at  CPA  is  then 
selected.  Sense  selection  ia  made  on 
the  first  logic  cycle  that  the  intruder 
ia  declared  a  threat.  Once  selected, 
the  sense  of  coaaand  for  a  specific 
intruder  does  not  change  throughout  the 
encounter. 

Since  the  •rose  is  selected  only  once, 
it  is  iaperative  that  the  occurrences  of 
incorrect  sense  choices  be  ttiniaised. 
To  determine  the  threat's  vertical  poai¬ 
tion  at  CPA,  the  BCAS  logic  projects  the 
threat's  vertical  position  at  tiae  of 
detection.  The  projection,  which  is 
based  on  the  tracked  vertical  rate  of 
the  intruder  at  time  of  detection,  may 
entail  as  much  as  a  35-second  position 
projection.  As  a  result,  che  peak  error 
in  che  cracked  vertical  race  at  tiae 
of  detection  can  be  multiplied  by  35  in 
che  determination  of  the  threat's 
projected  vortical  position  at  CPA. 
During  previous  testing  of  the  BCAS 
logic,  it  was  discovered  that  the  peak 
errors  in  the  vertical  race  estimate 
provided  by  the  a  -  I  cracker  often 
resulted  in  an  incorrect  sense  choice 
and  aarginal  separation  perforaance 
(reference  8). 
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The  third  use  of  the  tracked  vertical 
position  and  rate  is  aiade  in  the 
determination  of  the  severity  of  the 
BCAS  command.  The  logic  aelec'.a  one  of 
three  types  of  maneuvers:  positive 
commands  (climb  or  descend);  negative 
commands  (do  not  descend  or  do  not 
climb);  or  vertical  speed  limit  toam,ands 
(limit  descent  rate  or  limit  climb 
rate).  The  positive  commands  are  the 
most  severe,  since  they  require  a 
positive  response  by  the  BCAS  aircraft 
to  generate  separation.  Negative  and 
limit  commands  may  not  require  a  change 
in  the  BCAS  aircraft's  vertical  rate  in 
order  to  provide  proper  separation. 
Unlike  the  sense  of  the  command,  coeuaand 
severity  ran  tranaition  during  an 
encounter  poriod.  The  only  requirement 
is  that,  once  a  particular  command 
severity  is  selected,  it  most  be  dis¬ 
played  for  at  least  5  seconds  prior 
to  changing. 

Logic  testing  has  not  delected  a  command 
severity  selection  problem  chat  can  be 
traced  to  errors  in  vertical  rate 
estimation.  The  major  impact  on  large 
performance  caused  by  errors  in  vertical 
rate  estimation  occurs  in  the  area  of 
sense  selection. 
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PREVIOUS  VERTICAL  TRACKER. 

The  inputs  to  the  BCAS  vertical  tracking 
logic  are  the  mode  C  altitude  reports. 
These  reports  provide  altitude  infor¬ 
mation  to  the  nearent  100  feet.  The 
ideal  vertical  tracker  must  provide 
accurate  position  and  rate  estimates  in 
rpite  of  the  100-foot  granularity  in  the 
input  data.  Through  optimisation  of  the 
a  (altitude  tracking  constant),  and 
B  (rate  tracking  constant)  parameters, 
the  o  -  0  tracker  attempts  to  provide 
smooth  vertical  race  estimates  without 
excessive  delays  in  detecting  rate 
changes.  The  a  -  0  vertical  cracker 
used  by  Che  Active  BCAS  logic  originated 


with  the  intermittent  positive  control 
(1PC)  collision  avoidance  algorithms 
(reference  4 ) . 

With  consistent  1-second  updates,  thr 
vertical  poaition  estimate  of  thea-0 
tracker  at  tine  t  is  Z(t)  where 

Z(  t )  ■  ( 1- »)*(Z(  t-1 )  ♦  Z(t-l))  ♦a*7.r(t) 

(1) 

wh?re  Z(t)  and  Zf(t)  art  the  rate 
estimate  and  mode  C  report  at  tisui  t. 
The  altitude  tra..  ing  conatant  a  ■  0.4. 

The  vertical  rate  estimate  at  time  t  ia 

s 

Z(t)  and  ia  obtained  ea  follows: 

Z(t)  -  (l-0)*Z(t-l)  ♦  0*(Zr(t)  -  Z(t-l)) 

(2) 


where  the  original,  altitude  rate 
tracking  conatant  waa  9  ■  0.15. 

During  periods  of  missing  dats  reports 
(up  to  10  seconds)  (1)  is  replaced  by 

Z(t )  •  Z(t-At)  ♦  Z( t- At )*At  (3) 


where  the  duration  of  the  missing  data 
period  is  At. 

During  this  period,  the  rate  estimate  , - 
a  conatant,  Z(t-Ac).  The  initial  evalu¬ 
ation  of  the  BCAS  collision  avoidance 
logic  indicated  chat  the  probability  of 
incorrect  sense  choices  for  threats 
maneuvering  at  constant  vertical  rates 
waa  significant.  These  sense  choice 
problems  were  traced  to  the  errors  in 
Che  rate  estimate  at  the  time  of  threat 
detection  (references  5  and  6).  Since 
the  original  0  parameter  value  (0.15) 
was  based  on  the  4.7-second  updste  rate 
of  IPC,  the  presence  of  the  1-aecond 
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update  rate  of  Active  BCAS  cauaad  large 
quantitation  induced  arrort  in  the  rate 
eatiaate.  In  turn,  the  projection  of 
theae  large  rate  errora  tor  up  to  35 
aeconda  cauaad  an  incorrect  relative 
poaition  projection  at  CPA. 

Following  a  study  by  Broste  (reference 
7),  the  g  paraaeter  waa  reduced  from 
0.15  to  0.10.  Thia  change  ctuaed  the 
rate  eatiaate  to  be  aore  heavily 
aauothed  by  previoua  track  hiatory,  at 
the  coat  of  a  alight  loea  in  the 
tracker'a  tranaient  race  reaponae  tiae. 
Reference  7  deaonacratea  a  eignificant 
reduction  in  the  root  aear.  aquare  error 
in  the  rate  eatiaate  when  |  ia  reduced. 
However,  it  ia  iaportant  to  note  that 
aince  the  aenae  aelection  ia  the  result 
of  the  atochaatic  condition  of  the  rate 
eatiaate  at  a  aingle  inatant  in  tiaM,  it 
ia  the  control  of  the  maximal  error  in 
the  rate  eatiaate  that  ahould  be  of 
priaary  concern.  With  g»  0.10,  the 
cyclic  errora  in  the  rate  eatiaate  and 
che  pcaaible  error  in  the  vertical 
poaition  projection  lor  a  conatant  500 
feet  per  ainute  (fc/ain)  vertical 
aaneuver  are  ahown  in  figure  1.  The 
aode  C  report!  are  alao  ahown  in 
figure  1. 

At  a  conatant  500-ft/ain  vertical  rate, 
the  aode  C  report  will  change  once  every 
12  'cconda.  The  peak  quant i aat ion- 
induced  error  in  the  race  eatiaate 
exceeda  8  feet  per  aecond  (ft/aec). 
Thia  ia  the  aagnitude  of  the  vertical 
race  (500  ft/ain  ■  8.33  ft/aec).  The 
peak  error  in  the  race  eatiaate  occura  2 
aeconda  after  the  aode  C  report  changea. 
The  error  in  the  rate  eatiaate  reaulta 
in  a  cyclic  error  in  che  poaition 
projection.  During  che  12-second 
cracking  error  cycle,  che  error  in  the 
projected  poaition  exceeda  200  feet,  aix 
t iaea . 

Siailar  analyaia  of  che  cyclic  errora  in 
che  rate  eatiaate  cauaed  by  quantisation 
waa  performed  for  conatant  rate 
aaneuvera  between  500  and  3,500  fc/ain 
in  increawnta  of  100  ft/ain.  The  peak 


error  in  the  rate  eatiaate  and  the 
proportion  of  tiae  the  projected  poai¬ 
tion  error  exceeded  200  feet  are  ahown 
in  t able  1 . 

Table  1  indicatea  that  the  largeat 
percentagea  of  projected  poaition  error 
are  aaaociated  with  vertical  ratea 
between  500  and  1,100  ft/ain.  Nontero 
probabi litiea  of  the  projection  error 
exceeding  200  feet  appear  again  between 
1,600  and  1,900  ft/ain  and  between  2,700 
and  3,200  ft/ain.  Ratea  in  the  range  of 
1,600  to  1,900  ft/ain  arc  characterised 
by  tiaM  intervals  between  aode  C  transi¬ 
tions  which  (when  quantised  in  aeconda) 
do  not  result  in  conatant  time  measure- 
unnta  between  aode  C  report  transitions. 
At  i  .  500  ft/ain,  the  mode  C  report 
changes  every  4  aeconda;  it  changea 
every  3  aeconda  at  2,000  ft/ain.  There¬ 
fore,  ratea  between  1,500  and  2,000 
ft/ain  cause  mode  C  report  transitions 
to  cycle  between  3  and  4  aeconda.  It  ia 
thia  pattern  of  quantised  tiaM  aeasura- 
aenta  which  are  not  conatant  chat  cause 
the  peak  errora  in  che  rate  eatiaate. 
The  same  ia  true  in  the  2,700-  to  3,200- 
ft/ain  range.  Except  at  the  3,000- 
ft/ain  rate  where  the  transition  period 
is  a  conatant  2  aeconda,  each  of  che 
rates  represent  nonconatant  transition 
periods.  For  2,700-to  2 ,900~ ft/mi n 
rates,  the  transition  period  cycles 
between  2  and  3  seconds.  For  rates 
above  3,000  ft/ain,  it  cycles  between 
1  and  2  seconds. 

INITIAL  FLIGHT  TEST  RESULTS. 

Initial  live  flight  tests  were  conducted 
with  an  a  -  g  tracker  perforaing  the 
internal  BCAS  logic  vertical  tracking 
function.  The  paraaeter  values  were 
a  “  0.4  and  g  ■  0.1.  Alaost  immedi¬ 
ately,  probleaa  in  aense  choice  were 
detected  for  intruders  which  were 
essentially  in  level  flight.  The  sense 
choice  problems  were  traced  to  aode  C 
report  boundary  transitions.  The 
patterns  of  ax>de  C  reports  caused  by  the 
boundary  transitions  are  called  mode  C 
excursions . 
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TABLE  1.  QUANTIZATION- INDUCED  CYCLIC  PEAK  ERROR  IN  THE  RATE  ESTIMATE  AND  THE 

IMPACT  ON  THE  PROJECTED  POSITION  t  -  0.10 


RATE  PEAK  ERROR  PROPORTION  OF  TIME  PROJECTED 

(ft/min)  (ft/«ec)  POSITION  ERROR  EXCEEDS  200  FT 


500 

8.15 

0.50 

600 

7.96 

0.30 

700 

8.31 

0.23 

800 

8.03 

0.06 

900 

7.35 

0.09 

1,000 

5.78 

0.17 

1,100 

6.81 

0.09 

1,200 

4.90 

0.00 

1,300 

6.57 

0.07 

1,400 

4.67 

0.00 

1,500 

3.92 

0.00 

1,600 

5.72 

0.07 

1,700 

6.73 

0.05 

1,800 

4.90 

0.00 

1,900 

6.55 

0.09 

2,000 

2.84 

0.00 

2,100 

4.74 

0.00 

2,200 

5.62 

0.00 

2,300 

5.08 

0.00 

2,400 

3.19 

0.00 

2,500 

4.50 

0.00 

2,600 

3.82 

0.00 

2,700 

5.97 

0.06 

2,800 

6.41 

0.07 

2,900 

8.23 

0.07 

3,000 

1.61 

0.00 

3,100 

8.22 

0.07 

3,200 

6.40 

0.08 

3,300 

4.83 

0.00 

3,400 

6.21 

0.02 

3,500 

4.99 

0.00 
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Examples  of  mode  C  excursions  seen  in 
live  flight  testing  sre  shown  in 
figure  2. 


8 

ft. 


FIGURE  2.  LEVEL  FLIGHT  MODE  C  REPORT 
EXCURSION  PATTERNS 


In  the  first  case,  the  intruder  is  in 
level  flight  near  the  mode  C  transition 
boundary  of  3,050  feet.  A  slight 
deviation  in  the  flight  path  takes  the 
aircraft  above  the  transition  boundary 
causing  a  100-foot  change  in  the  quan¬ 
tized  altitude  report.  If  the  aircraft 
returns  to  below  3,050  feet  on  the  next 
cycle,  the  3,100-foot  mode  C  report  is 
not  reinforced.  On  the  other  hand,  the 
aircraft  can  remain  above  3,050  feet  and 
the  3,100-foot  mode  C  report  can  be 
reinforced.  The  quantization-induced 
rate  errors  for  these  two  mode  C 
excursions  are  shown  in  figure  3. 

Although  the  intruder's  flight  path  is 
nearly  level,  the  slight  altitude 
deviations  across  the  mode  C  transition 
boundary  can  induce  significant  errors 
in  the  rate  estimate.  For  the 
unreinforced  excursion,  the  peak  error 
of  10  ft/sec  occurs  when  the  mode  C 
report  changes.  The  rate  error  decays 
once  the  mode  C  report  returns  to 
its  previous  value.  The  peak  error  of 
10  ft/sec  could  cause  a  350-foot  error 


in  the  projected  vertical  position  In 
the  unreinforced  case,  the  error  in  the 
rate  estimate  is  depressed  below  the 
threshold  for  a  200-foot  error  in 
projected  position  1  second  after  the 
initial  excursion. 

For  cr sea  where  the  aircraft  remains 
above  the  transition  boundary  rein¬ 
forcing  the  mode  C  excursion,  the  peak 
error  in  the  rate  estimate  occurs  2 
seconds  after  the  mode  C  transition. 
The  peak  error  of  16.5  ft/sec  could 
cause  the  error  in  projected  position  to 
exceed  575  feet.  When  the  aircraft 
remains  abi  :e  the  transition  boundary, 
the  induced  errors  in  the  rate  estimate 
remain  above  the  200-foot  position  error 
threshold  for  10  seconds.  This  analysis 
reveals  that  with  a  -  8  tracking 
(«■  0.4  and  6*0.1)  mode  C 
quantization-induced  rate  errors  can 
remain  above  the  rate  threshold  that 
causes  200-foot  errors  in  position  pro¬ 
jection  for  1  to  10  seconds,  depending 
on  the  duration  of  the  mode  C  excursion. 

Both  mode  C  excursion  patterns 
investigated  occurred  often  in  flight 
testing.  Logic  changes  were  made  to 
control  the  impact  of  the  induced  rate 
error  caused  by  mode  C  excursions. 
Without  abandoning  a  -  0  cracking,  the 
following  CAS  logic  changes  were  made: 

1.  In  order  to  make  the  tracker 
less  sensitive  to  mode  C  transitions, 
the  rate  tracking  constant  0  was  reduced 
to  0.05. 


2.  The  sense  choice  logic  became  a 
"balanced"  logic,  in  the  sense  that 
neither  a  climb  nor  descent  avoidance 
maneuver  was  favored.  Previously, 
sense  choice  logic  had  favored  the 
selection  of  descent  sense  maneuvers. 

3.  A  rate  clip  parameter,  ZDLVL, 
was  added  to  the  sense  choice  logic.  If 
the  rate  estimate  was  less  than  10 
ft/second  any  sense  selection  that  might 
occur  on  that  logic  cycle  was  based 
strictly  on  relative  vertical  position. 
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FIGURE  3.  MODE  C  EXCURSION' INDUCED  TRACKED  RATE  ERRORS 


Subsequent  flight  testing  verified 
thst  these  logic  changes  eliminated 
sense  choice  problesis  associated  with 
■ode  C  excursions  for  aircraft 
essentially  in  level  flight.  However, 
the  reduction  of  6  to  0.05  has  signifi¬ 
cantly  reduced  the  tracker's  acceler¬ 
ation  detection  capabilities.  Broste 
(reference  7)  uses  the  term  "tracker 
equivalent  aesnry"  to  characterise  the 


tracker's  transient  rate  response  time 
as  a  function  of  a  and  0  .  Tracker 
equivalent  memory  approximates  the 
response  time  of  an  «  -  0  tracking 
system  to  input  transitions  and  identi¬ 
fies  the  tracker's  ability  to  detect 
accelerations.  Table  2  shows  the  cost 
incurred  in  terms  of  acceleration  detec¬ 
tion  with  the  reduction  of  0  to  0.05. 


TABLE  2.  TRACKER  RESPONSE  DELAY  AS  A  FUNCTION  OF  I 


Rate  Tracking  Constant 

_ | _ 


0.15 

0.10 

0.05 


Tracker  Response  Delay- 
Equivalent  Memory 
_ (Seconds) _ 


6.87 

9.06 

14.46 
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The  increase  in  smoothing  of  Che  rate 
estimate  with  the  reduction  of  0  to  0.05 
from  0.15  haa  resulted  in  the  doubling 
of  the  tracker’s  transient  rat?  response 
delay.  This  delay,  during  vertical 
accelerations  by  intruder  aircraft, 
could  critically  impact  BCAS  threat 
detection  by  inducing  delayed  alarms 
causing  insufficient  separation. 

NEW  NONLINEAR  TRACKING  APPROACH. 

In  order  to  provide  for  smoother  rate 
estimates  without  causing  a  signifi¬ 
cant  increase  in  the  tracker's  transient 
rate  response  delay  for  certain  condi¬ 
tions,  Lincoln  Laboratory  has  developed 
a  new  nonlinear  tracking  scheme.  The 
tracker,  which  is  considerably  more 
complex  than  the  a  -  0  tracker,  makes 
greater  use  of  the  mode  C  report  history 
than  the  a  -  0  tracker.  The  description 
of  the  tracker  is  presented  in  appendix 
A.  The  remaining  portions  of  this 
report  present  the  results  of  the 
comparison  of  the  nonlinear  tracking 
performance  with  the  a  -  0  tracking 
performance  (o  *  0.4,  0  ■  0.05). 


RESULTS  AND  ANALYSIS 


The  analysis  of  the  nonlinear  tracker 
performance  was  conducted  in  three 
phases.  The  fitst  characteristic 
analyzed  was  the  maximal  error  in  the 
rate  estimate  for  a  given  constant  rate 
vertical  maneuver.  To  avoid  the 
influence  of  tracker  initialisation 
procedures,  the  trackers  were  exercised 
with  the  fixed  vertical  rate  node  C 
inputs  until  the  truckers  reached  steady 
state  conditions.  The  resulting  cracker 
steady  state  conditions  may  contain 
oscillations  in  rate  estimation.  These 
oscillations  are  called  limit  cycle 
oscillations  (reference  9)  and  appear  in 
a  repeatable  pattern.  A  comparison  is 
made  between  the  maximal  errors  in  the 
limit  cycle  oscillations  for  each 
tracker.  In  this  phase,  mode  C  altitude 
reports  ere  provided  directly  to  the 


nonlinear  and  a  -  0  tracking  logic.  In 
the  second  phase,  each  tracker  was 
interfaced  into  the  Active  BCAS  logic; 
performance  was  analyzed  through 
simulation  using  the  Fast-Time  Encounter 
Generator.  In  the  final  phase,  selected 
live  flight  test  data  were  used  to 
compare  the  a  -  0  and  nonlinear  tracker's 
performance. 

STAND-ALONE  STEADY  STATE  PERFORMANCE. 

In  a  steady  state  500-ft/min  climb,  the 
mods  C  report  changes  every  12  seconds. 
Figure  4  presents  the  comparison  of  the 
limit  cycles  for  both  crackers.  This 
causes  the  periodicity  of  rate  error 
cycle  to  equal  12  seconds.  The  peak 
error  in  the  rate  estimate  within  that 
cycle  ave  the  errors  which  are  compared. 
Figure  5  present u  the  comparison  of  the 
peak  rate  error  for  both  trackers  for 
rates  ranging  from  100  to  3,500  ft/min. 
The  rate  error  which  would  yield  a  100- 
foot  vertical  position  projection  error 
(2.86  ft /sec)  is  also  identified. 

In  no  case  does  the  steady  state 
nonlinear  tracker  peak  rate  error  exceed 
the  a  -  0  tracker ' s  peak  race  error  for 
rates  less  than  400  ft/min.  Nominal 
errors  in  the  rate  estimates  which  occur 
for  both  trackers  are  large.  However, 
the  magnitude  of  the  rate  estimates  for 
both  trackers  is  less  than  10  ft/sec. 
Ac  a  result,  the  BCAS  logic  assumes  the 
aircraft  is  essentially  level  and 
ignores  the  rate  estimate  in  command 
sense  selection.  Current  relative 
position  is  used  instead  of  projected 
vertical  position  to  determine  command 
sense . 

The  altitude  bin  occupancy  time  is 
defined  as  the  duration  of  a  constant 
mode  C  report.  For  a  500-ft/min  rate, 
the  bin  occupancy  time  is  12  seconds. 
When  the  bin  occupancy  time  is  constant, 
the  nonlinear  tracker's  p*^ak  rate  error 
is  0  ft/sec.  Looking  at  figure  5,  this 
condition  occurs  for  rates  of  400,  500, 
600,  1,000,  1,200,  1,500,  2,000,  and 
3,000  ft/min.  This  occurs  because  these 
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FIGURE  4.  COMPARISON  OF  LIMIT  CYCLE  OSCILLATION  FOR  500  FT/M1M  CLIMB 


TllUt  VIATICAL  RATI  IPT/NM) 


FIGURE  5.  MAXIMAL  ERROR  COMPARISON  TOR  FINAL  STEADY  STATE  CLIMB 


races  represent  constant  inceger  bin 
occupancy  times  of  15,  12,  10,  6,  4,  3, 
and  2  seconds,  respectively. 

The  nonsero  peak  errors  in  the  nonlinear 
tracker  occur  for  rates  which  do  not 
yield  a  constant  bin  occupancy  time. 
For  instance,  a  2,700-ft/min  rate  1 45 
ft/sec)  yields  the  following  bin 
occupancy  sequence:  3  seconds, 
2  seconds,  2  seconds,  2  seconds,  3 
seconds.  The  nonconstant  bin  occupancy 
times  cause  the  rate  estimate  errors. 

Regardless  of  the  consistency  of  the  bin 
occupancy  time,  the  o  -  g  tracker  always 
has  a  nonsero  peak  rate  error.  Both 
trackers  have  mean  rate  errors  of  0 
ft/sec  for  all  steady  state  constant 
rate  conditions  analysed.  It  is 
important  to  note  that  the  constant  rate 
vertical  maneuvers  usually  employed  by 
pilots  (500,  1,000,  1,500,  and  2,000 
ft/min)  all  represent  constant  bin 


occupancy  times.  The  nonlinear  tracker, 
in  these  cases,  has  the  capability  of 
providing  error-free  rate  estimates. 

ACCELERATION  PERFORMANCE. 

The  stand-alone  analysis  was  extended  to 
characterise  tracker  performance  during 
and  following  vertical  accelerations. 
In  the  first  case  analysed,  mode  C 
inputs  represented  a  level  flight  air¬ 
craft  that  began  8  ft/sec2  (0.25  g) 
acceleration  at  time  t.  The  acceler¬ 
ation  continued  until  a  final  vertical 
rate  was  achieved.  The  analysis  was 
designed  to  compare  the  difference  in 
the  time  required  for  the  trackers  to 
obtain  an  accurate  estimate  of  the  final 
vertical  rate.  Conservative  estimates 
of  the  time  delay  were  obtained  by 
initially  having  the  level  flight  swde  C 
inputs  occur  while  the  aircraft  was  at 
the  lowest  portion  of  the  altitude  bin. 
In  this  way,  with  an  acceleration  into 
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the  climb  condition,  the  longest  period 
mi'ld  expire  prior  to  e  ■ode  C  report 
transition. 

In  figure  6,  the  final  rate  invaeti|ated 
wee  1,000  ft/min.  This  final  rate  is 
achieved  2.07  seconds  after  the  acceler¬ 
ation  begins.  The  first  aode  C  transi¬ 
tion  is  observed  at  t+7.  Rate  estimate 
for  the  a  -  I  tracker  is  6.48  ft/sec. 
The  nonlinesr  tracker  considers  this 
first  change  in  node  C  to  be  an  isolated 
change  since  the  previous  node  C  reports 
had  indicated  level  flight.  As  shown  in 
appendix  A,  the  nonlinear  rate  estimate 
at  t+7  is  8  ft/sec,  the  value  of  param¬ 
eter  PI.  The  next  mode  C  transition 
occurs  at  t+13.  B  i tween  t+7  and  t+13, 
the  a  -  0  rate  estimate  peaks  at  t+10. 
This  corresponds  to  the  occurrence  of 
the  peak  rate  estimate  for  reinforced 
■ode  C  excursion  rate  estimates  as  shown 
in  figure  3.  After  t+13,  the  nonlinear 
tracker  has  observed  the  constant  bin 
occupancy  time  for  the  1,000-ft/min  rate 
(6  seconds),  and  error-free  rate  esti¬ 
mates  occur.  The  a  -  0  rate  esti¬ 
mates  show  the  6-second  cyclic  pattern 
associated  with  the  bin  occupancy  tisie. 
The  e  -  0  rate  estimate  converges  to 
1,000  ft/min  but  will  oscillate  around 
that  value  indefinitely.  For  0.25  g 
acceleration  and  final  rate  of  1,000 
ft/min,  both  the  tracks  experience  the 
same  delay  in  providing  accurate  rste 
estimates  (13  seconds).  The  advantage 
of  the  nonlinear  tracker  is  the  error- 
free  condition  that  can  exist  after 
t+12. 

Figure  7  depicts  the  comparison  when  the 
final  rate  is  increased  to  2,000  ft/min. 
The  final  rate  is  achieved  in  4.15 
seconds.  The  first  mode  C  traneition 
occurs  at  t+6.  Aga  in,  this  initial 
transition  is  interpreted  as  an  isolated 
transition  by  the  nonlinear  tracker. 
The  seam  rate  estimate  (8  ft/sec)  occurs 
as  it  occurred  far  the  initial  transi¬ 
tion  in  the  1,000  ft/min  case.  The 
second  transition  occurs  at  t+9  and 
reflects  the  constant  3-second  bit 
occupancy  time.  For  times  greater  than 


t  +  9,  considerable  differences  in  the 
rate  estimates  provided  by  the  two 
trackers  exist.  The  cyclic  3-second 
pattern  of  rate  estimates  are  apparent 
in  the  a  -  0  rate  estimate.  Accurate 
estimates  of  the  rate  occur  10  seconds 
earlier  with  the  nonlinear  tracker.  The 
error  in  the  a  -  0  rate  estimate  exceeds 
10  percent  of  the  actual  rate  until 
t+19.  This  delay  in  obtaining  an 
accurate  estimate  of  the  final  rate  is 
due  to  the  large  equivalent  memory 
(table  2)  associated  with  the  0  ■  0.05 
value. 

In  figure  8,  results  for  a  final  rate  of 
3,000  ft/min  (2-second  bin  occupancy 
time)  are  shown.  In  this  case,  the 
final  rate  is  achieved  in  6.22  seconds. 
The  first  mode  C  transition  occurs  at 
t+5  and  the  second  at  t  +  7.  During  the 
acceleration  phase,  between  t  and  t+6, 
both  trackers  provide  very  poor  rate 
estimates.  However,  following  the 
second  transition  at  t+7,  the  nonlinear 
tracker  provides  an  error-free  rste 
estimate.  Due  to  the  equivalent  ammory 
of  the  a  -  0  tracker,  accurate  estimates 
of  the  true  rate,  50  ft/sec,  do  not 
occur  until  t+18.  Hence,  accurate  rate 
estimates  occur  11  seconds  sooner  with 
the  nonlinear  tracker. 

A  comparison  of  tracker  performance  for 
slow  accelerations  was  also  made.  The 
case  where  the  acceleration  rate  is 
0.1  g  (3.21  ft/sec2)  «nd  the  final 
rate  is  2,300  ft/min  will  be  reviewed. 
When  the  final  rate  is  2,300  ft/min,  the 
bin  occupancy  times  are  not  constant. 
The  bin  occupancy  pattern  is  a  combin¬ 
ation  of  2-  and  3-second  occupancy 
times.  The  results  of  the  compsrison 
are  shown  in  figure  9.  Achievement  of 
the  final  rate  required  11.92  seconds. 
Because  of  the  slow  acceleration,  the 
first  transition  does  not  occur  until 
t+8.  The  second  transition  does  not 
occur  until  t  +  11.  Although  the  non¬ 
linear  rate  error  exceeds  4  ft/sec 
between  t  +  11  and  t  +  15,  the  error  is 
lest  than  one-half  the  error  for  the 
a  -  0  tracker  for  the  same  time  period. 
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COMPARISON  OF  ACCELERATION  PERFORMANCE  —  ACCELERATION 
RATE  »  0.25  g,  FINAL  RATE  -  2,000  FT/MIN 
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FIGURE  8.  COMPARISON  OF  ACCELERATION  PERFORMANCE  —  ACCELERATION 
RATE  -  0.25  g,  Flt'AL  RATE  ■  3,000  FT/M1N 
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The  a  1 ov  acca lerat  ion  and  final  rata, 
which  did  not  reflect  constant  bin 
occupancy  times,  resulted  in  the  longest 
delays  for  tha  nonlinear  tracker  to  pro¬ 
vide  accurate  rate  estimates.  However, 
the  nonlinear  rate  estimates  were 
within  10  percent  of  the  true  rate  for 
tiaes  greater  than  L  ♦  1 5 .  Tha  same 
condition  did  not  result  for  tha 
•  -  I  tracker  until  t+25. 

Comparison  of  performance  during  high 
accelerations  (g  “  0.5)  was  also  made. 
Figure  10  presents  the  results  when  the 
final  rate  is  2,300  ft/ain.  With  the 
high  accelerat ion,  the  final  rate  is 
obtained  in  2.38  seconds.  The  first 
mode  C  transition  occurs  at  t+4. 
Although  the  rate  is  only  38.33  ft/sec, 
a  second  transition  occurs  at  t+6. 
Remember  that  the  2,300-ft/ain  rate 
reflects  a  combination  of  3-  and 
2-second  bin  occupancy  times.  Because 
the  second  transition  occurs  only  2 
seconds  after  the  first,  the  nonlinear 
rate  estimate  exceeds  the  true  rate  by  5 
ft/sec  between  c+6  and  ,.♦8.  The  magni¬ 
tude  of  the  error  is  leas  than  one-half 
the  error  in  the  a  -  6  rate  estiaete  for 
the  saae  period.  After  the  third  tran¬ 
sition  which  occurs  at  t+9,  the  errors 
in  the  nonlinear  rate  estimate  remain 
less  than  2.62  ft/sec.  Thus,  for  times 
greater  chan  t*9,  the  peak  error  in  Che 
nonlinear  tracker  is  leas  than  7  percent 
of  the  actual  rate.  This  same  accuracy 
does  noc  occur  in  the  a  -  0  rate  esti¬ 
mate  until  t+22.  Age  in,  the  nonlinear 
tracker  provided  accurate  rate  estimates 
12  seconds  sooner  than  t  h  e  a  -  ft 
tracker. 

The  preceding  analysis  reviewed  tracker 
performance  for  level  flight  aircraft 
which  accelerate  into  constant  rate 
climbs.  Siailar  results  can  be  expected 
for  transitions  frem  level  flight  into 
constant  rate  descents. 

The  comparison  of  the  performance  of  the 
nonlinear  and  a  -  0  trackers  during 
accelerations  from  level  flight  is  sum¬ 
marised  in  table  3.  If  the  acceleration 
is  initiated  at  time  t: 


1  t  +  * 

E*  “  —  ^|llvK)j  whara  *00  i*  the  rate 
k-t+i-10 
error  at  tiae  k. 

For  instance,  E10  is  the  average 
magnitude  of  the  tracked  rate  error 
during  the  first  10  seconds  following 
the  acceleration.  Similarly,  E^q  is  the 
average  magnitude  of  the  tracked  rate 
error  dur  ng  the  period  froa  10  tn  20 
seconds  after  the  initiation  of  the 
acceleration.  The  possible  everage 
magnitudes  of  the  vertical  position 
projection  error  can  be  obtained 
directly  by  multiplying  table  entries 
by  35  seconds.  The  coaparisons  are  made 
for  three  different  acceleration  rates 
and  final  vertical  rates  ranging  froa 
500  to  3,500  ft/ain. 

With  the  high  acceleration  rate 
(g  ■  0.5)  there  is  essentially  no 
difference  in  performance  during  the 
first  10  seconds  following  the  acceler¬ 
ation  when  the  final  rate  obtained  is 
1,500  ft/min  or  less.  However,  for 
final  rates  of  2,000  ft/min  and  abovj, 
the  nonlinear  average  rate  error 
approximates  one-half  to  two-thirds  of 
the  a  -  ft  tracker's  average  rate  error 
during  the  first  10  seconds  following 
the  acceleration. 

During  the  10-  to  20-second  time  period 
following  Che  acceleration,  the  non¬ 
linear  rate  errors  have  been  essentially 
eliminated  regardless  of  che  final  race 
obtained.  The  average  rate  errors  for 
the  a  -  ft  cracker  remained  significant 
during  Che  t+10  to  t+20  second  time 
period  when  final  rates  obtained  ware 
1,500  ft/min  or  more.  The  accuracy  of 
a  -  ft  rate  estimates  does  noc  match  che 
uniform  accuracy  of  the  nonlinear  race 
estimates  until  t*29  seconds. 

With  0.25  g  acceleration  rate,  the  saae 
pattern  of  results  is  observed.  For  the 
t  ♦  1 0  to  t+20  second  cime  period,  a 
larger  difference  in  performance 
results  than  was  observed  with  0.5  g 
accelerations.  Average  u  -  0  rate 
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TABLE  3.  COMPARISON  Of  TIKI  DEPENDENT  AVERAGE  ERROR  MAGNITUDES 
DURING  ACCELERATIONS  —  FT/SBC 


final  Rata 

Nonlinaar  Tracker 

a  - 

4  Tr  ack ar 

JL 

( ft/ain) 

hi 

hi 

hi 

*ia 

hi. 

0.50 

500 

8.33 

2.86 

1.65 

8.33 

2.31 

1.38 

1,000 

13.86 

2.34 

3.00 

13.51 

3.53 

1.29 

1,500 

16.41 

0.00 

C.00 

is.  03 

4  43 

0.99 

1 

2,000 

16.06 

0.00 

0 .  uC 

22.26 

5.04 

1.02 

2,500 

18.64 

1.43 

0.61 

27.26 

6.18 

1.40 

V 

3,000 

18.12 

0.00 

0.00 

32.34 

8,45 

1.28 

3,500 

18.09 

1.13 

0.77 

35.14 

9.48 

1.39 

0.25 

5*  0 

9.08 

3.43 

1.66 

9.08 

2.94 

2.64 

\ 

1,000 

13.75 

2.34 

o  on 

13.38 

3.52 

1.46 

1,500 

15.32 

0.00 

t'.CO 

17.71 

4.35 

0.90 

2,000 

13.66 

0.00 

Q.00 

21.70 

6.06 

1.02 

2,500 

18.09 

1.48 

0.67 

25.09 

7.46 

1.72 

3,000 

15.99 

0.00 

0.00 

28.31 

10.16 

1.51 

3,500 

18.42 

1.38 

0.65 

30.87 

13.19 

1.86 

0.10 

500 

7.63 

3.78 

2.15 

7.63 

2.:  3 

2.59 

1,000 

11.64 

4.43 

0.00 

11.88 

4.8V 

0.99 

1,500 

14.34 

1.92 

0.00 

14.46 

7.27 

1.54 

2  000 

15.52 

0.00 

0.00 

15.62 

10.38 

a  56 

1 

2,500 

15.52 

4.80 

0.95 

15.51 

15.19 

3.0s 

1 

3,000 

15.52 

8.17 

1.31 

15.51 

18.76 

4.92 

3,500 

15.52 

10.34 

1  30 

15.51 

22.16 

7.62 
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errors  for  the  t  +  10  to  t+20  second  time 
period  could  cause  vertical  projection 
errors  to  exceed  450  feet.  As  in 
the  high  acceleration  case,  the 
a  -  0  tracker  rate  accuracy  does  not 
approach  the  nonlinear  tracker  accuracy 
until  t+20  seconds. 

For  s low  accelerations  (g  ■  0.1),  the 
performance  of  Che  two  trackers  is 
nearly  identical  during  the  first  10 
seconds  following  Che  acceleration. 
This  result  occurs  independently  of  the 
final  rate  obtained.  When  the  final 
rate  obtained  exceeded  1,000  ft/min,  the 
nonlinear  tracker  performed  considerably 
better  than  the  a  ~  8  tracker  during  the 
10~  to  20-second  period  following  the 
acceleration.  For  slow  accelerations, 
the  o  -  8  tracker  rate  errors  persist 
during  the  t  +  20  to  t  +  30  second  time 
period . 

With  a  -  8  tracking,  projection  errors 
as  large  <.s  265  reet  could  still  persist 
more  than  20  seconds  after  the 
acceleration. 

DECELERATIONS  TO  LEVEL  FLIGHT 
CONDITIONS. 

In  this  section,  we  will  review  tracker 
performance  when  an  aircraft  in  a  fixed 
rate  climb  returns  to  level  flight.  The 
climb  rates  investigated  ranged  from  500 
to  3,500  ft/min.  The  analysis  of  the 
tracker  performance  during  the  return  to 
level  flight  from  a  climb  rate  of  2,300 
ft/min  will  be  presented.  This  rate 
represented  the  severest  conditions 
analyzed . 

In  figure  11,  a  -0.1  g  deceleration 
effect  is  presented.  At  -0.1  g,  it 
takes  11.92  seconds  to  return  to  level 
flight  from  a  2,300-ft/min  climb  rate. 
The  mode  C  transition  pattern  of  2-  and 
3-second  bin  occupancy  times  has  been 
established  in  the  climb.  The  detection 
of  the  deceleration  requires  the  tracker 
to  detect  a  deviation  from  the  2-  and 
3-second  bin  occupancy  pattern.  At  t+6, 
the  nonlinear  tracker  detects  tne 
excessive  bin  occupancy  time  following 


the  deceleration.  As  a  result,  the  rate 
estimate  is  adjusted  with  a  decaying 
process.  The  same  procedure  is  used  at 
t+9.  The  adjustment  almost  eliminates 
the  error  in  the  nonlinear  rate  esti¬ 
mate.  The  final  mode  C  transition 
occurs  at  t+6.  At  t  +  14,  the  nonlinear 
tracker  observes  that  the  estimated  bin 
occupancy  time,  ZM0D7,  has  been  exceeded 
by  5  seconds  (the  value  of  parameter 
P5).  As  a  result,  the  rate  estimate  is 
set  to  0  ft/sec  t  u  represent  level 
flight.  The  a  -  8  rate  error  exceeds 
5.71  ft/sec  for  an  additional  e  seconds. 
During  this  period,  the  a  -  8  errors 
would  result  in  at  least  a  20)-foot 
error  in  vertical  position  projection. 

For  a  higher  deceleration  rate  (figure 
12),  the  performance  of  the  nonlinear 
tracker  is  even  better  when  compared  to 
the  a  -  8  tracker.  For  -0.5  g  deceler¬ 
ation,  an  aircraft  climbing  at  2,300 
ft/min  can  level  off  in  2.38  seconds. 
The  last  mode  C  transition,  prior  to 
deceleration,  occurs  at  t-1.  Henc >,  at 
t+3  the  nonlinear  tracker  detec  s  an 
excessive  bin  occupancy  time  and  adjusts 
the  rate  estimate.  At  t+6,  the  esti¬ 
mated  bin  occupancy  time  his  been 
exceeded  by  5  seconds  and  the  rate 
estimate  is  reset  to  0  ft/sec.  The 
a  -  8  rate  error  does  not  decay  below 
5.71  ft/sec  until  t+12.  Again,  we  see  a 
6-second  difference  in  the  establishment 
of  accurate  rate  estimates. 

During  decelerations  to  level  flight, 
the  nonlinear  tracker  resets  the  rate 
estimate  to  0  ft/sec  when  the  excess  bin 
occupancy  time  exceeds  5  seconds.  The 
a  -  8  tracker  cannot  reset  its  rate 
estimate  to  0  ft/sec.  The  8  value  of 
0.05  causes  the  decay  in  the  rate  esti¬ 
mate  to  0  ft/sec  to  be  quite  slow. 
Depending  on  the  precision  of  the  logic 
arithmetic  operations,  the  a  -  8  rate 
estimate  may  cycle  indefinitely  with  a 
nonzero  magnitude.  This  could  cause 
problems  with  the  CAS  logic's  vertical 
speed  limit  routines.  These  routines 
make  explicit  checks  for  nonzero 
vertical  rates. 
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FIGURE  11.  DECELERATION  TO  LEVEL  FLIGHT  CONDITION  PERFORMANCE  g  -  -0.1 


[•j 


Table  4  summaiizen  performance  of  che 
two  trackers  during  and  following 
decelerations  from  established  climbs  to 
level  flight.  The  comparison  is  made 
for  three  deceleration  rates  and  for 
established  climb  rates  ranging  from  500 
to  3,500  ft/min.  For  almost  all  condi¬ 
tions  tested,  the  nonlinear  tracker 
declared  level  flight  within  10  seconds 
of  the  initiation  of  the  decelerations. 
The  large  errors  that  still  exist  in  the 
nonlinear  rate  estimate  more  than  10 
seconds  after  initiation  of  the  deceler¬ 
ation  occurred  for  the  high  initial 
climb  rate  (3,000  and  3,500  ft/min)  and 
for  the  slow  deceleration  rate  (-0.1  g). 
In  these  cases,  level  flight  conditions 
do  not  occur  until  15.5  and  18.1  seconds 
after  the  deceleration  was  initiated. 
The  a  -  B  rate  estimates  contained  a 
considerably  larger  error  than  the  non¬ 
linear  tracker  during  the  10-  to 
20-second  time  period  following  the 
initiation  of  deceleration. 

DECELERATION  TO  OTHER  THAN  LEVEL  FLIGHT 
CONDITIONS . 

Analysis  of  tracker  performance  during 
decelerations  from  a  fixed  rate  to 
another  rate  (not  level  flight)  was 
made.  In  figure  13,  the  initial  rate  is 
-2,500  ft/min  and  at  time  c  a  deceler¬ 
ation  to  -1,000  ft/min  is  initiated. 
With  a  -0.25  g  deceleration,  the  new 
rate  is  obtained  in  3.10  seconds.  At 
t+3,  the  nonlinear  tracker  detects  an 
excess  in  the  bin  occupancy  time  and 
adjusts  the  rate  estimate.  Between  t+3 
and  t  +  5,  the  nonlinear  tracker  over¬ 
estimates  the  deceleration.  However, 
the  magnitude  of  the  rate  error  of  the 
nonlinear  track  is  less  than  the  a-B 
tracker  error  magnitude.  Beyond  t+6, 
the  nonlinear  tracker  rate  estimate  is 
error  free.  Consistently  accurate  a-B 
rate  estimates  do  not  occur  until  t+14, 
an  8-second  difference  in  the  delay. 

PITCHOVER  MANEUVERS. 

Maneuvers  in  which  a  constant  rate 
deceleration  takes  an  aircraft  from  a 


constant  climb  rate  into  a  constant 
descent  rate  were  analysed.  In  figure 
14,  the  aircraft  is  established  in  a 
3,000-ft/min  climb  rate,  then  deceler¬ 
ates  into  a  1,000-ft/min  descent.  The 
deceleration  rate  presented  is  -0.5  g. 
The  nonlinear  tracker  detects  excess  bin 
occupancy  time  at  t+3  and  adjusts  the 
rate  estimate  to  approximately  1,000 
ft/min.  At  t+7,  a  downward  transition 
in  altitude  bins  is  detected  and  a 
negative  rate  estimate  occurs.  The 
a  -  B  tracker  does  not  provide  a  nega¬ 
tive  rate  estimate  until  t+13.  At  this 
time,  the  nonlinear  rate  estimate  is 
error  free.  In  general,  the  analysis  of 
pitchover  maneuver  results  indicated 
significantly  earlier  (6  to  10  seconds) 
negati*  c  rate  detection  by  the  nonlinear 
tracker.  The  differences  in  time,  to 
obtain  accurate  final  rate  estimates, 
approximated  differences  identified  in 
the  comparison  of  deceleration 
performance. 

During  pitchover  maneuvers  the  ability 
of  the  vertical  tracker  to  determine 
when  the  direction  of  vertical  movement 
has  changed  is  critical.  The  parameter 
ZDLVL  is  used  to  discount  low  vertical 
races  by  intruder  aircraft  in  selecting 
maneuver  sense.  The  delay  in  che 
tracker  detecting  e  rate  less  than  ZDLVL 
is  a  good  measure  of  tracker  performance 
during  pitchover  maneuverr.  Tracker 
responses  during  end  following  -0.5  g 
pitchover  maneuvers  were  analyzed. 
Initial  climb  races  varied  between  500 
and  3,500  ft/min.  The  -0.5  g  pitchover 
was  sustained  until  a  1,000-ft/min 
descent  rate  was  established. 

Table  5  identifies  the  results  of  the 
pitchover  analysis.  The  comparison  is 
made  in  terms  of  the  delay  in  seconds 
until  certain  tracked  conditions  exist. 
Since  the  delay  periods  are  dependent  on 
the  location  within  an  altitude  bin  when 
the  pitchover  maneuver  is  initiated,  the 
delay  is  stated  as  a  period  of  time. 
Since  the  maximum  rate  estimates  for 
both  trackers  are  less  than  10  ft/sec 
for  a  500-ft/min  climb,  the  tracked 
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TABU  4. 


wwr™1WH  ur  TINS -DEPENDENT  AVERAC1 
WRING  DECELERATIONS  TO  LEVEL  FLICH1 


ERROR  MAGNITUDE 


Initial 
Cl i*b  Rate 
8  (ft/min) 


“0.50 

500 

1,000 

1.500 

2,000 

2.500 
3,000 

3.500 

“0.25 

5  00 

1,000 

1.500 

2,000 

2.500 
3,000 

3.500 

-Q.  10 

500 

1,000 

1.500 

2,000 

2.500 
3,000 

3.500 


Nonlinear  Tracker 

ho  ho 


6.76  0.0  0.0 

13.67  0.0  0.0 

7.61  0.0  o.O 

12.51  0.0  0,0 

5.52  0.0  o.O 

9.40  0.0  0.0 

18.74  0.0  0.0 


6.46  0.0  0.0 

12.08  0.0  0.0 

5.63  0.0  0.0 

8.60  0.0  0.0 

12.72  0.0  0.0 

16.16  0.0  0.0 

20.10  0.0  0.0 


5.79  0.0  0.0 
10.68  0,48  0.0 
11. 24  0.61  0.0 

11.71  0.70  0.0 
13.62  0.51  0  0 
10.83  6. 71  0.0 

10.71  9.90  0.18 


a-  •  Tracker 

*10  *20  ho 


4.10  0.64  0.08 
13.84  2.89  0.45 
12.65  2.34  0.30 
19.87  4.11  0.46 
20.86  4.36  0.58 
25.31  5.78  0.76 
32.41  8.83  1.20 


4.07  0.70  0.08 
12.98  3.25  0.45 
10.76  2.34  0.30 
16.45  4.21  0.55 
22.39  7.03  0.96 
23.90  8.93  i.25 
28.07  12.12  1.72 


3.40  0.74  0.07 
10.39  3.25  0.45 
11.24  5.49  0.80 
11*71  7.92  0.98 
13.62  10.32  2.06 
10.83  16.08  5.00 

10.71  19.81  6.94 
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FIGURE  14.  P1TCH0VER  MANEUVER  PERFORMANCE  COMPARISON 
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TABLE  5.  COMPARISON  Or  TRACKER  RESPONSE  DELAYS  FOLLOWING 
-0.5  g  PITCHOVER  MANEUVER  TO  -1,000  FT/M1N 


Period  Tracked  Rate  Period  Until  Descent 


Initial 

Climb  Rate 
(ft/ain) 

Exceeds  10  ft/sec 
(seconds) 

Nonlinear  a-S 

is  Detected 
(seconds) 
Nonlinear 

o-S 

500 

- 

- 

4-6 

5-7 

1,000 

2-7 

5-9 

5-8 

8-13 

.  ^0 

3-6 

6-10 

6-8 

8-13 

2,000 

3-6 

7-10 

4-7 

10-14 

2,500 

4-7 

8-11 

6-8 

12-14 

2,700 

4-7 

8-11 

6-8 

12-14 

3,000 

4-8 

9-11 

6-8 

13-15 

3,500 

5-8 

10-12 

6-8 

13-16 
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race*  do  not  exceed  10  ft/aec  during  the 
picchover  when  the  initial  cliab  rate  ia 
500  ft/min.  The  aaxiaua  tiae  required 
to  depraae  the  nonlinear  rate  eatiaate 
below  10  ft/sec  following  the  initiation 
of  the  pitchover  ia,  in  general,  leaa 
chan  the  ainiaua  delay  chat  ia  expected 
with  a  -  ft  tracking.  The  delay  with 
«  -  0  tracking  required  before  a  deacent 
is  detected  ia  Cwice  the  delay  experi¬ 
enced  with  nonlinear  cracking.  Thia 
reault  ia  fairly  uniform  acroaa  all 
inicial  cliab  races  above  \  ,000  ft/min. 

STAIRSTEP  MANEUVER. 

The  final  acand-alone  tracker 
coapariaona  analysed  atairatep  vertical 
profilea.  The  aircraft  waa  assumed  to 
be  deacending  at  \  conatant  rate.  A 
level-off  maneuver  waa  then 
accomplished,  followed  by  deceleration 
into  a  descent  5  to  20  seconds  later. 
Figure  15  presents  a  typical  vertical 
profile. 


The  results  are  shown  in  figure  16.  The 
SMjor  advantage  of  the  nonlinear  tracker 
is  the  ability  to  obtain  an  accurate 
level  flight  rate  estimate  much  sooner 
tnan  the  o  ••  6  tracker.  The  nonlinear 
tracker  provides  the  level  flight  rate 
eatimatea  for  the  time  period  between 
t+10  and  t+25.  The  n  -  6  tracker 
obtains  accurate  level  flight  rate  eati¬ 
matea  only  for  times  greater  than  t*2Q. 
Once  the  deacent  ia  resumed  at  t+26, 
the  nonlinear  tracker  responds  as  if 
the  aircraft  had  always  been  in  level 
flighc.  As  long  aa  the  level  flight 
period  exceeds  the  altitude  strata 
occupancy  time  (3  seconds  for  a 
2 , 000“ f t /mi n  vertical  rate)  by  5 
seconds,  the  nonlinear  tracker  will 
reset  the  vertical  rate  estimate  to 
0  ft/min  for  the  level  flight  portion  of 
the  stairstep  maneuver. 


FIGURE  15.  TYPICAL  STAIRSTEP  MANEUVER 
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FIGURE  16.  STAIRSTEP  MANEUVER  PERFORMANCE  COMPARISON 


FAST-TIME  ENCOUNTER  GENERATOR  RESULTS. 

Previous  problems  with  command  sense 
choice  during  intruder  icce ieret ions 
were  identified  in  reference  8. 
Unequipped  intruder  encounter  scenarios 
like  the  sequence  depicted  in  figure  17 
resulted  in  BCAS  commands  which  reduced 
vertical  separation  at  CPA.  Fast-time 
aiaulation  analysis  was  conducted  10 
see  ;.f  the  earlier  detection  of  acceler¬ 
ations  by  the  nonlinear  tracker  had 
iaproved  the  collision  avoidance  logic 
pertoraance  for  scenarios  siailar  to 
chose  shown  in  figure  17. 

The  variety  of  geoaetries  analysed  are 
shown  in  figure  lb.  The  Active  BCAS 
logic  does  not  perait  sense  of  coseaand 
to  change  once  it  is  selected.  If  the 
wrong  sense  is  selected  and  the  logic 
selects  a  positive  conaand,  vertical 
separation  at  CPA  is  reduced  by  the 
incorrect  sense  choice.  The  response 
model  used  for  this  analysis  included  a 
5-second  pilot  response  cieliy . 

The  question  that  is  addressed  in  this 
section  is  "How  long  does  the  error  in 
the  vertical  estimate,  which  wo»>U  cause 
the  wrong  command  sense  to  be  selected, 
persist  following  the  level  off  maneuver 
by  the  intruder?"  For  conditions  shown 
in  figure  17.  the  incorrect  sense  choice 
for  the  BCAS  aircraft  is  climb.  Since 
sense  choice  can  occur  on  any  logic 
cycle  when  all  threat  detection  condi¬ 
tions  are  satisfied,  the  improvement  in 
sense  choice  because  of  nonlinear 
vertical  tracking  can  be  identified  by 
comparing  the  period  of  incorrect  sense 
choice  associated  with  each  cracker 
following  the  level  off  maneuver. 

Table  6  presents  the  reduction  in 
duration  of  the  incorrect  sense  choice 
period  which  occurred  with  nonlinear 
cracking.  Following  the  level-off 
maneuver,  both  the  a  -  8  and  nonlinear 
tracking  rate  estimates  contain  a 
residual  error.  For  a  specific  plsnned 
separstion  value,  an  incorrect  sense 
choice  results  when  the  residual  rate 


error  projected  over  35  seconds  exceeds 
the  planned  vertical  separation  value. 
The  improvement  in  sense  choice  per¬ 
formance  with  nonlinear  tracking  is  due 
to  a  faster  tracker1*  response  to  the 
variation  in  the  mode  C  report  pattern. 

Table  b  indicates  that  sense  choice 
improvement  increases  with  the  magnitude 
of  the  initial  vertical  rate.  For 
larger  values  of  vertical  separation, 
larger  residual  errors  sre  necessary  for 
incorrect  sense  choices  to  result.  For 
lower  initial  vertical  rates  (less  than 
-1400  ft/min)  and  the  larger  planned 
vertical  separr.t  ions ,  the  duration  of 
incorrect  sense  choice  periods  for  both 
trackers  are  nearly  identical.  However, 
for  higher  initial  vertical  rates, 
nonlinear  tracking  reduced  the  incorrect 
sense  periods  by  up  to  7  seconds. 

Table  7  presents  the  improvement  with 
nonlinear  cracking  that  results  when 
the  level-off  maneuver  is  performed  more 
rapidly  (-0.5  g  versus  -0.25  g). 

Although  a  considerable  reduction  in  the 
duration  of  the  incorrect  sense  choice 
period  occurs  with  nonlinear  tracking, 
wrong  sense  choices  may  still  occur. 
Wrong  sense  choices  result  with  the 
nonlinear  tracker  because  .he  sense 
selection  logic  does  not  make  optimum 
use  of  the  nonlinear  tracker.  Research 
should  be  conducted  to  improve  the  use 
of  nonlinear  tracking  data  by  the  sense 
choice  logic.  Two  immediate  consider¬ 
ations  should  be  investigated.  Since 
command  sense  cannoc  change  once  it  is 
selected,  when  sufficient  time 
(-range/range  rate  sufficiently  large) 
until  CPA  remains,  command  sense  selec¬ 
tion  should  be  delayed  during  periods  of 
acceleration  by  the  intruder  Several 
internal  variables  in  the  nonlinear 
tracker  (ZM0D8  or  a  function  of  ZM0D7) 
appear  to  be  good  indicators  of  acceler¬ 
ation  by  the  intruder.  A  second,  more 
sophisticated  approach  is  to  base  sense 
determination  on  vertical  position 
projection  which  includes  a  vertical 
acceleration  factor. 
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FIGURE  18.  GEOMETRIES  USED  IN  COMPARISON  OF  DETECTION  PERFORMANCE 


REVIEW  OF  LIVE  FLIGHT  TEST  RESULTS. 

During  Active  BCAS  flight  testing,  en 
encounter  occurred  with  e  DC-9  aircraft. 
The  Federal  Aviation  Adainietrat ion 
(FAA)  BCAS  aircraft,  registration  number 
N49,  waa  level  at  16,000  feet.  The  aode 
C  report  history  revealed  that  the  DC-9 
was  cliabing  at  approxiaate ly  2,700 
ft/ain.  About  35  seconds  prior  to  tha 
tracked  ainiaua  range  tiae,  the  sade  C 
report  history  showed  a  vertical 
deceleration  by  the  DC-9.  The  highest 
aode  C  report  observed  was  15,500  feet. 
The  deceleration  rate  was  about 
(-0.5  g).  Tha  deceleration  period 
lasted  approximately  4  seconds.  The 
aode  C  report  history  shows  a  final 
vertical  rate  of  about  -1,000  ft/ain. 


after  the  aode  C  history  shows  a 
deceleration  by  the  DC-9.  At  this 
tiae,  the  projected  vertical  separation, 
Zp,  is  calculated  as  follows: 

Zp  ■  BCAS  altitude  -  (DC-9  tracked 
position  ♦  35  seconds  x  DC-9 
tracked  altitude  rare) 

For  the  o  -  I  tracker  in  the  airborne 
experimental  BCAS  unit 

Zp  -  16,000  -  (15,543  ♦  35*(34)) 

■  -593  ft 

This  implies  the  DC-9  is  projected  to  be 
593  feet  above  the  level  flight  BCAS 
aircraft  at  CPA.  As  a  result,  the  BCAS 
logic  displayed  a  "do  not  cliab"  cosaiand 
for  the  BCAS  aircraft. 


Tiie  vertical  profile  for  this  encounter 
is  shown  in  figure  19.  The  initial  BCAS 
alara,  do  not  cliab,  occurs  4  seconds 
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This  encounter  was  simulated  using  the 
fast-tiae  encounter  generator.  The  tiae 
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TAAL8  6.  RIDUCTION  IN  DURATION  OF  INCORIICT  8IXN8B  CHOIC*  PIRIOD 
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nATA  RECORDED  AND  PROCESSED 
BY  THE  FAA  TECHNICAL  CENTER 


TABLE  7.  REDUCTION  IN  DURATION  OF  INCORRECT  SENSE  CHOICE 

PERIODS  IN  SECONDS  ACHIEVED  WITH  NONLINEAR  VERTICAL 
TRACKINC  (-0.5  g  DECELERATION) 


Initial 
Vertical 
RATE  (Ft/Min) 
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FIGURE  19.  VERTICAL  PROFILE  OF  DC-9  ENCOUNTER 


versus  altitude  plot  of  the  results 
which  occurred  with  a  -  0  tracking  is 
shown  in  figure  20.  In  simulation,  the 
sense  of  the  initial  alarm,  which  was 
descend,  was  incorrect.  In  simulation, 
command  selection  appears  to  have 
occurred  1  second  later  than  during  the 
actual  flight.  If  a  first  hit  had 
occurred  at  the  48th  second,  the 
command  selected  would  have  been  "do 
not  climb"  since  the  absolute  valuO  of 
VMD  |-544|  >470  feet,  the  threshold  for 
positive  commands.  The  logic  requires  a 
command  to  be  displayed  for  a  minimum  of 
5  seconds  prior  to  being  replaced  by 
another  comswnd.  A  review  of  the  live 
flight  data  indicated  a  positive 
descent  command  did  not  result  because 
the  projected  vertical  separation  was 
less  than  470  feet  during  the  5-second 
period  immediately  following  the  display 
of  the  initial  "do  not  climb"  alarm. 
The  a  -  8  tracker  rate  estimates  for  the 
DC-9  remained  positive  for  12  seconds 
following  the  vertical  deceleration. 
During  this  period,  the  collision 


avoidance  logic  would  sense  the  positive 
late  as  vertical  closure  and  possibly 
select  an  incorrect  command  sense. 

By  repeating  the  fast-time  simulation, 
results  were  obtained  for  nonlinear 
tracking.  Figure  21  shows  the  results. 
With  the  nonlinear  tracker,  a  negative 
vertical  rate  estimate  occurs  5  seconds 
after  the  deceleration.  Since  sense 
selection  occurs  at  this  time,  the  sense 
selected  climb  is  correct.  Since  the 
projected  vertical  miss  distance  of  670 
feet  exceeds  the  positive  command 
threshold  of  470  feet,  a  "do  not 
descend"  command  is  displayed. 

The  tracker  responses  associated  with 
this  encounter  were  investigated  further 
by  supplying  the  DC-9~mode  C  history 
directly  to  the  trackers  independent  of 
the  collision  avoidance  logic.  The 
tracker  responses  are  shown  in 
figure  22.  Assuming  the  deceleration 
wao  initiated  at  time  t,  the  duration 
of  the  periods  of  incorrect  sense  choice 
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DATA  recorded  and  processed 

BY  THE  FAA  TECHNICAL  CENTER 
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FIGURE  20.  LIVE  FLIGHT  TEST  ENCOUNTER  —  ALPHA-BETA  TRACKING  RESULTS 


DATA  RE  CORO 
BY  THE  FAA  TE 


FIGURE  21.  LIVE  FLIGHT  TEST  ENCOUNTER  —  NONLIRER  TRACKING  RESULTS 


VERTICAL  RATE  (FT/ SEC) 


l- 


F1GURE  22.  TRACKER  RESPONSES  FOR  DC-9  VERTICAL  PROFILE 
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for  each  cracker  war*  obtained.  The 
aenae  choice  logic  haa  a  vertical  rate 
threshold,  ZDLVL,  auch  that  aenae  choice 
ia  baaed  atrictly  on  current  relative 
position  when  the  ablolute  value  of  the 
race  estimate  ia  leaa  than  ZDLVL 
(10  ft/sec).  For  the  DC-9  vertical 
profile,  the  duration  of  the  incorrect 
aenae  choice  period  equals  the  time  it 
cakea  the  tracker  to  depress  the  rate 
estimate  below  10  ft/sec.  With  the 
nonlinear  tracker,  thia  occurs  at  t+4; 
hence,  the  incorrect  sense  choice  period 
is  3  seconds.  The  a  -  8  tracker  does 
not  depress  the  rate  estimate  below 
10  ft/aec  until  t+9,  resulting  in  an 
8-second  incorrect  sense  choice  period. 
This  corresponds  closely  with  the  live 
flight  test  results.  The  rate  estimate 
in  the  flight  teat  results  remained 
above  10  ft/sec  for  10  seconds  following 
the  DC-9  deceleration. 

U'l  the  nonlinear  tracker,  the 
stand-alone  analysis  indicated  a  level 
flight  condition  was  declared  6  seconds 
after  the  deceleration  was  initiated. 
The  e  -  8  tracker  required  12  seconds  to 
establish  the  sme  result.  This  dif¬ 
ference  indicates  the  a  -  8  tracker 
would  declare  a  vertical  closure  condi¬ 
tion  for  6  more  seconds  than  the 
nonlinear  tracker. 


CONCLUSIONS 


1.  Review  of  a  -  8  cracker  performance 
through  the  use  of  faat-time  simulation 
results  in  the  following  conclusions: 

a.  The  previous  rate  tracking 
constant  value  of  8a  0.15  had  been 
developed  for  use  with  4.7-second  data 
rates.  The  much  faster  data  rate  (1 
second)  of  Active  Beacon  Collision 
Avoidance  System  (BCAS)  caused  large 
errors  to  be  induced  into  the  rate 
estimate  by  mode  C  quantisation.  The 
large  errors  often  resulted  in  incorrect 
command  sense  choice  even  during  steady 
state  vertical  maneuvers  by  the  threat 
aircraft. 


b.  A  reduction  in  8  to  0.1 
somewhat  reduced  the  occurrence  of 
incorrect  sense  choice  during  steady 
state  vertical  maneuvers.  However, 
initial  flight  testing  showed  that  100- 
foot  mode  C  excursions  for  aircraft  in 
level  flight  could  result  in  rata  errors 
exceeding  16  feet  per  second  (ft/sec). 

c.  Modifications  to  the  sense 
choice  logic  and  a  reduction  in  8  to 
0.05  eliminated  sense  choice  problems 
associated  with  100-foot  mode  C  report 
excursions.  However,  the  decrease  in 
8  has  resulted  in  a  significant  increase 
in  the  tracker's  transient  rate  response 
delay . 

2.  Review  of  the  nonlinear  tracker 
described  in  the  appendix  led  to  the 
following  conclusions: 

a.  The  complexity  increase  when 
compared  to  the  e  -  8  tracker  requires 
the  storage  of  ten  additional  logic 
parameters.  The  number  of  elements  in 
the  own  aircraft  state  vector  increased 
by  eight.  Also,  the  number  of  elements 
in  each  intruder  threat  file  increased 
by  eight. 

b.  The  missing  or  garbed  altitude 
report  flag  in  the  Collision  Avoidance 
System  (CAS)  logic  ZFLG  was  easily 
incorporated  into  the  nonlinear  tracking 
logic  to  indicate  when  missing  altitude 
report  logic  should  be  used. 

3.  Comparative  analysis  of  the 
performance  of  the  nonlinear  tracker 
with  the  •  -  8  tracker  resulted  in  the 
following  conclusions: 

a.  For  100-foot  mode  C  excursions, 
the  maximal  rate  error  associated  with 
nonlinear  tracking  was  8  ft/sec.  This 
occurs  because  transitions  in  the  mode  C 
reports  for  aircraft  with  tracked 
vertical  rates  of  0  ft/sec  are  treated 
as  isolated  transitions.  The  maximal 
error  occurs  at  the  time  of  transition 
and  decays  even  when  successive  reports 
reinforce  the  transition.  The  maximal 
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•rror  in  the  rate  eatiaate  ia  leaa  then 
ZDLVL.  Hence,  incorrect  aenae  choicea 
do  not  occur  dur  in*  100-foot  mode  C 
report  excursions  with  nonlinear 
tracking. 

b.  During  ateady  atate  condition! 
for  all  ratea  analysed,  the  nonlinear 
tracker  peak  error  in  the  rate  estimate 
was  aignificantly  aaaller  than  the  a-  | 
tracker  peak  error  in  the  rate  eatiaate. 

c.  During  ateady  atate  conditions 

for  vertical  ratea  which  repreaent 
conatant  bin  occupancy  tines,  the  non¬ 
linear  tracker  would  provide  error-free 
rate  eatiaates.  The  rates  associated 
with  constant  bin  occupancy  tiaes  (500, 
1,000,  1,500,  2,000  feet  per  minute 

(ft/ain))  are  ratea  at  which  aircraft 
noraally  maneuver.  Regardless  of  rate, 
the  a  -  0  tracking  rate  estimate  cycles 
with  a  periodicity  equal  to  the  bin 
occupancy  tinea. 

d.  Following  acceleration  from 
level  flight  into  conatant  rate  cliaba, 
little  difference  in  the  tiaes  required 
for  each  tracker  to  obtain  accurate  rate 
estimates  is  detected  for  final  vertical 
rates  £1,000  ft/min.  However,  an  the 
final  rate  ia  increased,  nonlinear 
tracking  results  in  accurati  rate 
estimates  significantly  earlier  (6  to  12 
seconds)  than  with  the  a  -  0  tracker. 
This  fact  appears  uniform  across  the 
acceleration  rates  tested. 

e.  The  weakei.t  performance  by  the 
nonlinear  tracker  occurs  during  slow 
accelerations  to  final  rates  which  are 
not  associated  witn  constant  bin  occu¬ 
pancy  times.  However,  accurate  race 
estimates  are  obtained  considerably 
earlier  with  the  nonlinear  tracker  than 
with  the  «  -  0  tracker. 


f.  During  and  following  level-off 
maneuvers,  the  nonlinear  tracker 
required  significantly  leas  time  (6  to 
10  seconds)  to  obtain  0  ft/aec  rate 
estimates.  When  the  excess  bin  occu¬ 
pancy  time  exceeds  5  seconds,  the  rate 
estimate  ia  reset  to  0  ft/sec.  The 
a-  0  tracker's  rate  estimate  can 
oscillate  around  0  ft/sec  inuefinitely , 
depending  on  the  arithmetic  precision  of 
the  logic. 

g.  During  pitchover  maneuvers 
earlier  detection  of  the  change  in  the 
sign  of  the  vertical  rate  occurs  with 
the  nonlinear  tracker.  The  large 
transient  response  delay  associated  with 
0  “  0.05  results  in  long  periods  in 
which  the  e  -  0  tracking  rate  estimate 
has  the  wrong  sign. 

h.  During  stairstep  maneuvers 

whenever  the  level  flight  duration 
exceeds  5  seconds  plus  the  expected  bin 
occupancy  time,  the  nonlinear  rate  esti¬ 
mate  is  reset  to  the  correct  value  of 
0  ft/sec.  This  does  not  occur  with 
*  -  0  tracking. 

i.  Fast-time  simulation  of  Active 
RCAS  logic  with  both  nonlinear  anda-0 
tracking  has  shown  that  a  significant 
improvement  in  command  sense  choice 
occurs  with  nonlinear  tracking.  Addi¬ 
tional  improvement  in  sense  choice 
performance  during  periods  of  acceler¬ 
ation  may  be  obtained.  This  improvement 
may  be  possible  by  augmenting  the  sense 
selection  logic  with  logic  chat  accounts 
for  accelerating  intruders. 

4.  Review  of  live  flight  test  sense 
choice  anomalies  indicates  that  the 
occurrence  of  incorrect  st.r'se  choice 
would  be  reduced  with  nonlinear 
cracking. 
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RECOMMENDATIONS 


Based  on  the  results  of  the  analysis 
of  the  nonlinear  tracker,  it  is 
recommended  thst  the  review  tracker  be 
incorporated  into  Active  Beacon 
Collision  Avoidance  System  (BCAS)  flight 
testing.  Furthermore,  it  is 

recoHaended  that  additional  analysis  be 
made  to  identify  methods  of  optimising 
Collision  Avoidance  System  (CAS)  logic 
sense  selection.  The  optimisation 
should  make  use  of  the  additional 
information  that  is  available  with 
nonlinear  cracking  that  was  not 
previously  available  with  Alpha-Beta 
Tracking. 
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APPENDIX  A 

NONLINEAR  VERTICAL  TRACKER 


Only  minor  changes  were  made  to  the  rode 
provided  by  Lincoln  Laboratory  in  the 
VTRK66  Version  of  the  nonlinear  altitude 
tracker.  The  changes  included:  (1)  the 
parameter  P2  was  deleted  since  it  was 
not  referenced  in  the  code,  (2)  the 
previous  a  -  I  tracking  of  own  altitude 
was  replaced  with  the  nonlinear  tracking 
logic,  and  (3)  the  check  for  missing  or 
unreliable  altitude  reports  (ZMa0)  was 
changed  to  conform  with  the  flag  pro- 
vided  by  Beacon  Collision  Avoidance 
System  (BCAS)  surveillance  tracking 
which  indicates  missed,  garbled,  or  an 
unreliable  altitude  report,  ZFLAGa0. 
This  change  permitted  the  direct  inter* 
face  of  VTRK66  into  the  Collision 
Avoidance  System  (CAS)  logic  intruder 
tracking  routiue,  TR1ACT. 

To  support  the  addition  of  VTRK66, 
specific  parameters  had  to  be  added  to 
the  BCAS  CAP  logic  parameter  list.  The 
parameters  added  are  shown  in  table  A-l. 

Two  other  parameters  uaed  in  VTRK66 
logic  were  not  needed  by  the  CAS  logic. 
The  parameter  DT,  which  identifies  the 
nominal  update  rate  (1  second)  already 
exists  within  the  CAS  tracking  modules. 
The  BCAS  variables  TDATA  and  TDATAI 
already  determine  the  duration  of 
missing  report  periods.  The  quanti¬ 
sation  bin  with  parameter  Q  was  not 
added  to  the  parameter  list.  The 
constant  100  was  added  in  its  place  in 
the  VTRK66  logic.  Two  parameters  were 
deleted  from  the  BCAS  parameter  list. 
They  were  the  o  -  8  position  tracking 
constant,  ALFAZ,  and  the  rate  tracking 
constant,  BETAZ. 

VTRK66  logic  required  that  several  new 
elements  be  added  to  the  intruder  track 
file.  Tne  elements  are  shown  in 
table  A-2. 


Each  element  in  table  A-2  ia  subscripted 
with  the  intruder's  ID  or  track  file 
number.  The  tame  elements  were  also 
added  to  the  own  aircraft  state 
vector.  It  is  important  to  note  that 
2MOD6  was  added  to  the  intruder  state 
vector  although  the  same  information 
already  exists  in  the  state  vector  in 
the  form  of  TDATA.  This  was  done 
because  ZM0D6  is  internally  manipulated 
in  the  VTRK66  logic  and  may  have  an 
effect  on  other  usea  of  T)ATA  in  the  CAS 
logic . 

Certain  additional  track  file  elements 
now  require  initialization  prior  to  use 
by  the  tracker.  When  the  BCAS  ia  first 
made  operational,  own  intruder  state 
elements  must  be  initialized.  The  addi¬ 
tional  elements  requiring  initialization 
are  shown  in  table  A-3. 

When  intruder  reports  are  first  passed 
to  the  CAS  logic  by  the  surveillance 
logic,  procedures  must  initialize 
certain  new  elements  in  the  intruder 
track  file.  The  initialization  is  per¬ 
formed  in  TR1ACT.  Nonlinear  tracking 
occurs  on  subsequent  reports.  The  ele¬ 
ments  of  the  intruder  state  vector  to 
be  initialized  are  shown  in  table  A-4. 

When  an  intruder  track  file  exists  and 
surveillance  tracking  provides  at 
least  a  range  report  on  a  particular 
logic  cycle,  the  nonlinear  tracking 
logic  as  shown  in  figure  A-l  is  used  to 
provide  intruder  tracking.  This  logic 
will  handle  unreliable,  missing,  or 
garbled  altitude  reports.  When  no 
surveillance  report  ia  received  for  an 
intruder  with  an  established  track  file, 
coasting  of  the  track  with  the  non¬ 
linear  cracker  is  accomplished  as  shown 
in  figure  A-2.  Own  cracking  is 
accomplished  by  integrating  the  non¬ 
linear  tracker  into  TKOACT  in  the  same 
manner  shown  in  figure  A-l.  The  only 
difference  is  that  the  missing  altitude 
report  logic  can  be  deleted. 


TABLE  A-l.  BCAS  CAS  LOGIC  PARAMETER  LIST  ADDITIONS 


Parameter 

Definition 

Nominal  Value 

PI 

Magnitude  of  rate  allowed  following  isolated 
altitude  transition 

8  feet  per  second 

P3 

Rate  decay  factor  when  nc  reinforcing 
transition  has  occurred 

0.90 

PA 

Stiff  rate  saoothing  parameter 

0.0A 

P5 

Excess  bin  occupancy  time  which  results  in 
transition  to  level  flight 

5  seconds 

P6 

Excess  bin  occupancy  time  which  results  in 
correction  to  rate  estimate 

1.5  seconds 

P7 

Amount  of  discrepancy  in  bin  occupancy  time 
which  triggers  rate  reinitialization 

1.5  seconds 

P8 

Parameter  used  to  position  an  estimated  bin 
transition  within  an  interval  of  missing  reports 

0.6 

P9 

Position  smoothing  parameter 

0.3 

P10 

Decay  factor  for  residual  monitor,  ZMODIO, 
for  scans  without  detection  of  excess  residual 

0.0 

Pll 

Bin  occupancy  smoothing  parameter  when  excess 
residuals  are  detected 

0.6 

P12 

Decay  factor  following  a  correction  to  ZMODIO 

0.75 

P13 

Value  of  ZM0D9  at  which  transition  from  start 
up  smoothing  to  normal  smoothing  occurs 

18 
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TABLE  A-2 . 


ELEMENTS  ADDED  TO  INTRUDER  TRACK  FILE 


VTRK66 

Element 

Previous  Track 
File  Element 

Definition 

Units 

,ZM 

ZR1NT 

Mode  C  report.  ZFLAGO  when  report 
is  missing  or  garbled 

feet  (100) 

ZMOD1 

ZINT 

Tracked  altitude  of  intruder 

feet 

ZMOD2 

ZD  I  NT 

Tracked  altitude  rate  of  intruder 

feet  per  second 

ZMOD3 

New  Element 

Time  last  mode  C  report  was 
received  for  this  intruder 

seconds 

ZMOD4 

New  Element 

Previous  mode  C  report  for 
this  intruder 

feet  (100) 

ZMOD5 

New  Element 

Time  of  transistion  to  previous 
mode  C  altitude  for  this  intruder 

seconds 

ZMOD6 

New  Element 

Time  of  last  CAS  track  update 
for  this  intruder 

seconds 

ZMOD7 

New  Element 

Estimated  bin  occupancy  time 

seconds 

ZMOD8 

New  Element 

Firmness  of  rate,  typically  equal 
to  number  of  altitude  transitions 
observed  for  current  rate 

ZMOD9 

New  Element 

Start  up  counter  needed  to 
establish  rate 

— 

Z MOD 10 

New  Element 

Residual  indicator.  Used  to  detect 
a  trend  in  tracker  residuals  which 
indicate  vertical  acceleration 

— 
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TABLE  A' 3. 


OWN  STATE  VECTOR  INITIALIZATION  PROCEDURE 


VTRK66  Element 

Previous  CAS  Logic  Element 

Initialization  Value 

ZMODl 

ZOWN 

Curient  own  mode  C  report 

ZM0D2 

ZDOWN 

0  feet  per  second** 

ZMOD3 

New  Element 

TCUR* 

ZNODA 

New  Element 

Current  mode  C  report 

ZM0D5 

New  Element 

Added  below 

ZM0D6 

New  Element 

TCUR 

ZM0D7 

New  Element 

100/ (ABS (ZMOD2 )+0 . 1 ) 

ZN0D8 

New  Element 

5 

ZN0D9 

New  Element 

0 

ZMODl 0 

New  Element 

0 

ZM0D5 

New  Element 

TCUR-ZMOD7  ♦  7 

ZM 

ZKOWN 

Current  mode  C  report 

♦Variable  name  for  current  BCAS  time  in  CAS  logic. 

**If  BCAS  is  turned  on  in  flight  there  should  be  an  initial  standby  period. 
This  will  permit  the  possible  observance  of  a  bin  occupancy  period.  The 
standby  period  should  last  5  seconds  or  until  two  mode  C  transitions  are 
observed,  whichever  is  shorter.  If  less  than  two  transitions  occurred, 
ZM0D2  should  be  initialized  to  0  ft/sec.  If  two  transitions  are  observed, 
then  ZM0D2"(mode  C  report  on  second  transition  —  mode  C  report  on  first 
transition)/tine  between  transitions.  This  procedure  will  reduce  the  time 

delay  required  to  obtain  accurate  own  aircraft  rate  estimates  when  own 
aircraft  is  maneuvering  at  high  vertical  rates  (above  20  ft/sec)  when  BCAS 
is  first  initialized. 
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TABLE  A-4.  INTRUDER  STATE  VECTOR  INITIALIZATION  PROCEDURES 


VRTK66  Element 

Previous  CAS  Track  File  Element 

Initialisation  Value 

ZMODl 

ZINT 

Surveillance  altitude 
report  ZRINT 

ZM0D2 

ZD1NT 

Surveillance  rate  ZDSV 

ZN0D3 

New  El«Mnt 

TCUR 

ZN0D4 

New  Element 

ZRINT 

ZM0D5 

New  Element 

Added  below 

ZM0D6 

New  Element 

TCUR 

ZN0D7 

New  Element 

100/(ABS(ZDSV)-0. 1) 

ZH0D8 

New  Element 

5 

ZH0D9 

New  Element 

0 

ZMODl 0 

New  Element 

0 

ZMODS 

New  Element 

TCUR-ZM0D7  ♦  7 

ZM 

ZRINT 

ZRINT 
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OT*TCUR-TDATAI 

DTN«TCUR-Z* 

ZP-Z1*DTM*Z2 


DZN* ZM-Z4 
DBINS*ABS(DZH)/100 
Z*=Z**1  -M0»D«IMS 


Zf LAC*  0 


HISSING  ALTITUDE  REPORT  LOGIC 


Z1  =  ZP 


z*sph 


START-UP  SMOOTHING 

21  =  ZP*  Pf*(ZM-ZP) 
P1ZJ»Z2*P4* 

(ZM-ZP)/DTM 


DZM-  0 


DZH«  0 


|  TRANSISTION  LOGIC 


Z4=  ZM 
ZS-TCUR 


ZS  -TCUR*FI*  L.HO 
(Z3-TCUR-M)  r 


Z3>  Zt 


B  J  NO  TRANSISTION 
LOGIC  / 


ZFLAG =1 


Z3-TCUR 


ZA-TCUR 


CONTINUE  WITH 
INTRUDER  RANGE 
TRACKING 
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FIGURE  A-l.  NONLINEAR  TRACKING  LOGIC  Aa  INTEGRATED  INTO  TRIACT  (Sheet  l  of  3) 
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FIGURE  A-l. 


NONLINEAR  TRACKING  LOGIC  AS  INTEGRATED  INTO  TRIACT  (Sheet  2  of  3) 
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SURVEILLANCE  TRACKER  DECLARES  MISSING  REPORT 


DELETE  INTRUDER 
STATE  VECTOR 


DT=TCUR- 
ZP-Z1-*-  (T 
Z1  =  ZP 
Z6=TCUR 

TDATA 

CUR-Z6)  *Z2 

R  -  R  +  DT*RD 

TDATA  =  TCUR 

COAST  ALTITUDE 


COAST  RANGE 
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FIGURE  A-2 .  INTEGRATION  OF  NONLINEAR  TRACKING  LOCIC 
INTO  MISSING  REPORT  LOGIC  OF  TRIACT 
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